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Abstract 

The structure of the electric field in a plasma has been elucidated for the skin 
effect problem. An expression for the distribution function in the half-space and the 
electric field profile have been obtained in the explicit form. The absolute value, the 
real part, and the imaginary part of the electric filed have been analyzed in the case 
of the anomalous skin effect near to a plasma resonance. It has been demonstrated 
that the electric field in the skin effect problem is predominantly determined by 
the discrete spectrum, i.e., the oscillation frequency of external field is the value of 
plasma frequency. 

1. Introduction. Statement of problem. 

The skin effect is associated with the response of an electron gas (in a 
metal or in a gas plasma) to an external alternating electromagnetic field 
that is tangential to the surface P, El- This classical problem has been 
studied by many authors ||3j - [B] and, up to now, has remained the subject 
of investigation. The main attention has been focused on the calculation 



of the impedance. The distribution function of electrons and the electric 
field in plasma almost have not been investigated previously. 

It has been demonstrated that the electric field is the sum of the integral 
term and two (or one) exponentially decreasing particular solutions to the 
initial system and that one particular solution disappears depending on 
the anomaly parameter. 

Let's Maxwell plasma fills the half-space x > 0. Here x is the orthogonal 
coordinate to the plasma boundary. Let's the external electric field has 
only y component. Then the self-consistent electric field inside in plasma 
also has only y component Ey(x,t) = E(x)e~^'^*. We consider the kinetic 
equation for the electron distribution function: 

| + .,| + eE(.)e-g = .(/o-/(t,..v)). (1) 

In (1) z/ is the frequency of electron collisions with ions, e is the charge 
of electron, /o(v) is the equilibrium Maxwell distribution function, p = mv 
is the momentum of electron. 



/o(v)=nr^ exp(-/?V), p = 



Try 2kBT 
Here m is the mass of electron, ks is the Boltzmann constant, T is the 
temperature of plasma, v is the modulus of the velocity of the electron, n 
is the concentration of electrons (number density), c is the speed of light. 



The electric field E{x) satisfies to the equation: 

E"{x) = / Vyfit, X, v) d^v. (2) 

We assume that intensity of an electric field is such that linear approx- 
imation is valid. Then distribution function can be presented in the form: 

/ = /o (1 + C'y exp{-iujt)h{x, 12)) , 

where C = ^/^y is the dimensionless velocity of electron, ji = Cx- Let 
/ = vtt is the mean free path of electrons, vt = vt is the thermal 

electron velocity, r = l/iy. We introduce the dimensionless parameters 
and the electric field: 

= yt, Xi = y, e(xi) = 

Later we substitute xi for x. The substitution produces the follow- 
ing form of the kinetic equation (1) and the equation on a field with the 
displacement current (2): 

dh 

'^dx ^ ^ ^^^^ 

00 

e"(x) = —i-^ / exp(—iJ.'^)h(x,iJ.')diJ.'. (4) 
V TT J 

—00 

Here, 



100 ^ A 

a = -rT-, zq = 1 — \l = ujT = —, = 



2 



6 is the classical depth of the skin layer III, ctq = , ctq is the electric 

conductance, a is the anomaly parameter. 

The boundary conditions at the plasma surface for the distribution func- 
tion of electrons in the case of specular reflection of electrons from the 
boundary can be written as follows [T]: 

/i(0, /i) = /i(0, -/i), < /i < oo. (5) 

The distribution function will be sought in the form of a decaying func- 
tion far from the boundary; that is, 

/i(+oo,/i) = 0, -oo < /i < 0. (6) 

The electric field deep in the plasma far from the surface decays. Taking 
into account this circumstance, the boundary conditions for the electric 
field are written in the form 

e(0) = 1, (7) 

e(+oo) = 0. (8) 
2. Decomposition on eigenfunctions 

The separation of variables in (3) and (4) within several steps leads to 
the exponentially decreasing solutions 

hrj{x, Id) = exp(-^)$(T/, Id), erj{x) = exp{-^)E{ri), (9) 



where the separation parameter (also termed the spectral parameter) r] 
continuously fills the interval (0, oo), which, therefore, is called the contin- 
uous spectrum of the problem. 

Substitution of relationships (9) into the initial system of equations (3) 
and (4) leads to the characteristic system of equations 

(r/-^)$(r/,/i) = -E{r]), 

00 

Oi [ 

T J 

— oo 

The functions ^{rj^ii) and £"(77), which are referred to as the eigen- 
functions of the characteristic system and correspond to the eigenvalue (or 
characteristic value) of the parameter 77, are defined by the expressions 

a 1 

$(r/, /i) = ^rfP h Xif]) exp(r]^)6(r] - ji), (10a) 

Vtt 7/ - /i 

£(„) = (106) 



where the dispersion function \{z) (see, for example, [3]) is given by the 
formula 

+00 

^, , (^z'^ f exp(— /i^) , ia ,^ , 

X(z) = l + ^ / — ^d/i, a= ^. (11) 

Vtt J /i - z 

—00 

With the use of the argument principle, it is possible to show that, in 
the {a, Q) plane, there exists a domain (Fig. la) so that, if the point 
a,n) G D^, the dispersion function has four zeros ±770 and ±771, and 



if {a,Q) G D~ (where D~ is exterior of the domain D~^), the dispersion 
function has two zeros ±770- The designations rjQ and 771 correspond to 
the zeros with the positive real parts: Re 770 > and Re 771 > 0. The 
boundary of the domain Z)+ is found from the equation a;^(/i) = and, in 
the parametric form, is determined by the equations 1 — 30,^ ± a q{ii) = 

0, sri - ri^ - ap{id) = 0, -00 < ii< +00. 

It should be noted that the parameters a and Q, are proportional to the 
electric field frequency; i.e., they are not independent. In this respect, it 
seems quite natural to introduce the dimensionless independent frequencies 



Ui = , l/i = 



UpVc LJp Vc 



and to construct the corresponding domains and (Fig- lb) in their 
plane. Here, a = a;i/(z/f), O = uji/iyi, cUp is the plasma frequency and n 
is the electron concentration. 

The zeros 770 and 771 correspond to the following eigenfunctions of the 
characteristic equation that are associated with the discrete spectrum: 

Hm, ^) = V ^(^^) = ^ = 0' 1- 

The zeros of the dispersion function can be calculated in the explicit 
form with the use of the formulas for its factorization. In the case of two 



zeros ±?7o, the dispersion function (see [3j) can be represented in the form 



\{z) = a{rf,-z^)X{z)X{-z), 



where 



X{z) =expy(z), 

oo 

1 [\liG{T)dT 





In the case of four zeros ±770 and ±771, the dispersion function can be 
written as follows: 

\{z) = a(r/2 - z^M - z^)X,{z)X,{-z), 

where the function Xi[z) is expressed through the function X[z): Xi[z) = 
X{z)/{z-l). 

By calculating the left- and right-hand sides of the former formula for 
the factorization of the dispersion function (for, example, at the point 
2; = 0), after some transformations, we obtain the relationship for its zeros 

1 

±?yo = = + 1. 

^aX{z)X{-z) 

In the skin effect theory, the normal and anomalous skin effects are 
recognized [7]. In the case of the normal skin effect, the mean free path 
of electrons is considerably smaller than the skin depth; i.e., the anomaly 
parameter satisfies the inequality a <C 1. The anomalous skin effect cor- 



responds to the case where the mean free path of electrons is considerably 
larger than the characteristic skin depth: a ^ 1. 

Let us construct the general solution to the initial system of equations 
in the form of the expansion in eigenfunctions of the discrete and con- 
tinuous spectra. Since the discrete spectra for zero and unit indices are 
different and the continuous spectrum does not depend on the index, the 
expansions of the solution in both cases differ only in the nonintegral terms 
corresponding to the discrete spectrum. 

In [3], it was demonstrated that the distribution function of electrons 

and the electric field, which are the solution to the problem described by 

expressions (3)-(8), have the following expansions: 

1 



7/ N a y;^ AkTjl f ZqX 

h{x, ii) = —=y exp 



+ / exp ( ) A{T^)^{T^,ii)dT], (12) 







e{x) = Akfii exp ( - ^ ) + 

fc=o ' 



oo 



"""^^ ^ exp(-^)r/2^(r/)dr/. (13) 



"n J \ T] 



Here, Re r]k > 0, Ak {k = 0, 1) are unknown constant coefficients of 
expansions (12) and (13) (the so-called coefficients of the discrete spec- 



trum), and A{r]) is an unknown function (the so-called coefficient of the 
continuous spectrum). 

It should be noted that, in the case of two zeros of the dispersion func- 
tion, it is necessary to set A; = in relationships (12) and (13). Therefore, 
the structure of the electric field depends on the domain that contains 
the point with the parameters (a, fl). 

In [H|, it was shown that the coefficient A(r]) of the continuous spectrum 
is represented in the form 



zo/A+(r/)A-(7y) 
where 

oo 



27T J X{ir) 

—oo 

The coefficients of the discrete spectrum are written in the following 
form: 

^ , A: = 0,1. 

azoIr]lX'{r]k)' 

The impedance is given by the formula 

^ 4:7Tiu e(0) 

According to the boundary conditions for the filed, we have e(0) = 1. 
Therefore, the following expression holds true for the impedance: 



4:7Tiuj STTicol 

" ^V(0) " 



oo 



1 f dr 



7T J X{ir) 







3. Distribution function and the electric field 

With the use of the determined coefficients of the continuous and dis- 
crete spectra, the electric field profile in the half-space can be represented 
in the explicit form 

1 / ZoX\ 1 / ZoX\ 

oo 


Formula (14) will be subsequently used for analyzing the behavior of 
the electric field in the half-space. 

Now, we consider the profile of the distribution function of electrons in 
the half-space in the explicit form. The distribution function is represented 
in the form of two terms: 

h{x, ji) = hd{x, ji) + hc{x, ji), 

where the terms hd{x^ ji) and hc{x, fi) correspond to the discrete and contin- 
uous spectra, respectively. With the use of the equality for the coefficients 
of the discrete and continuous spectra, these terms are written as follows: 



I, I \ 1 / ( «oa:\ »;exp(-»)2) 
= exp ( - 





At the plasma boundary, i.e., at x = 0, the last relationship can be 
calculated in the explicit form. As a result, we have 



hc{0,id) 



1 



1 



1 



27rizoI J VA+(?7) [ij) J ij — jj, 





X{fi)6{jj.) exp(;U-^ 



1 



2izQlay/Tiy?' \X'^{fi) X {fi) 



V^), /i G (-oo,+oo). 



Here, 6{ii) = latO<;U<oo and 6{ii) = at — oo < n < 0. 



By using contour integration methods [5j , for the first term we find that 



/ic(0,/i) 



ZqI 



+ 



^0 



/iA(/i)(l-^(/i)) 

X+{fi)X~{fi) ' (r/o - /i)A'(r/o) 



1 

+ - 



r^dr 



{iji - fl)X'{r]i) 7T J X{iT){T'^ + 11^ 



By summing up the terms corresponding to the discrete and continuous 



spectra in the distribution function, we finally obtain 



/i(0,/i) = 



ZqI 



oo 



dr 



X+{fi)X-{fi) ttJ X{it){t^ + fi^] 





From this expression, for the distribution function of electrons moving 
at the metal boundary (i.e., in the case — oo < < 0), we have 



ZqI 



00 



1 

+ - 



dr 



A+(/i)A-(/i) 7tJ X{it){t^ + 11^) 





1 < /i < 0, 



For electrons specularly reflected from the metal boundary, we derive 



00 



/i(0,/i) = 



r^dr 



ttzqI J X{iT){T'^ + 





-, < /i < 1. 



These functions satisfy the specular boundary condition h{0, fi) = h{0, — /i). 
For all subsequent figures, we consider the typical case with the ratio 
^^/c = 0.003. 

The behavior of the real and imaginary parts of the distribution function 
at the boundary is illustrated in Fig. 2. In view of the specular boundary 
condition, the distribution functions of electrons reflected from the bound- 
ary (0 < /i < oo) and electrons moving to the boundary (— oo < /i < 0) are 
symmetric with respect to the point /i = 0. The functions are constructed 
for the parameters o; = 1, Q = 333. Figure 2a depicts the real part of the 
electric field, fig. 2b depicts the imaginary part of the electric field. Let us 
compare the imaginary part of the electric field Im ec{x) with the real part 
of the field Re ec{x). The electric-field amplitude | Im (/i(0, //))| > 6 • 10^ 
is considerably bigger for the imaginary part. It should be noted that the 
distribution function rapidly decreases with an increase in the quantity /i. 
This circumstance is a manifestation of the ineffectiveness concept [S], ac- 
cording to which only electrons moving almost parallel to the surface, i.e., 
for which the quantity /i is considerably smaller than unity, are significant 
in the case of the anomalous skin effect. 

The real and imaginary parts of the electric field in the vicinity of the 
boundary are presented in Figs. 3. The curves depicted in Fig. 3 cor- 



respond to the following parameters: e: = 10 ^ (a = 900, Q, = 1000) for 
curves 1,£ = 3 • 10"^ (a = 100, n = 333) curves 2, and £ = 9 • 10"^ (a = 11, 
Q, = 111). This is an anomalous case. All the curves are considered near 
plasma resonance, i.e. the value 7 = 1 and u = Up. 

Figure 3a shows the real part of the electric field Re ed{x), which corre- 
sponds to the discrete spectrum. An increase in the anomaly parameters 
leads to a drastic decrease in the depth of penetration of the electric field 
deep into the electron plasma. 

Figure 3b depicts the real part of the electric field Re ec(x), which cor- 
responds to the continuous spectrum. In this case, the electric-field am- 
plitude — I Re ec(x)| < 1, 2 • 10~^ is considerably smaller that that for the 
real part due to the discrete spectrum. 

The real part of the electric field Re ec(x), which is associated with the 
continuous spectrum, is nine orders of magnitude smaller than the real 
part of the field Re ed{x) corresponding to the discrete spectrum. 

Therefore, the real part of the electric field in the vicinity of the plasma 
boundary is actually determined by the discrete spectrum. 

The imaginary part of the electric field Im ed{x), which is associated 
with the discrete spectrum, is shown in Fig. 3c. As the anomaly parameter 
increases, the depth of penetration of the imaginary part of the electric 



field deep into tlie plasma decreases slowly in contrast to the depth of 
penetration of the real part. 

Figure 3d presents the imaginary part of the electric field Im ec{x) 
which corresponds to the continuous spectrum. Let us compare the imag- 
inary part of the electric field Im ec{x) with the real part of the field 
Re ec{x). It can be seen that the amplitude has the same order of magni- 
tude: I Im ec(x)| < 1, 2 • 10~^. Therefore, the imaginary part of the electric 
field at the aforementioned values of the parameter, in actual fact, is also 
determined by the discrete spectrum. However, the imaginary part of the 
electric field corresponding to the discrete spectrum is nine orders of mag- 
nitude bigger than the real part of the electric field corresponding to the 
discrete spectrum. 

As can be seen from the plots presented in Fig. 3, the contribution of 
the discrete spectrum at the aforementioned values of the parameter to 
the electric field is considerably larger than the contribution of the contin- 
uous spectrum. Thus, the above analysis of the electric field strength has 
demonstrated that, in the case of the anomalous skin effect, the electric 
field strength is determined in the vicinity of the boundary by the discrete 
spectrum. 

Figure 4 the modulus of the electric field in the case of the anomalous 



skin effect. On the X axis is taken logarithmic scale. The curves depicted in 
these figures correspond to the following parameters: e = 3- 10~^ and 7 = 5 
{a = 5,n = 1666) for curve 1, e = 3 • 10"^ and 7 = 1 (a = 100, n = 333) 
for curve 2. It can be seen from Fig. 4 that with anomaly parameter 
increases in 20 times the modulus of the electric field in the anomalous 
case decreases one order of magnitude more rapidly. 
Conclusions 

Thus, in this paper, we have demonstrated that the electric field and the 
distribution function of electrons for the skin-effect problem are determined 
by their particular solutions. These solutions are the sums of the solutions 
corresponding to the discrete spectrum (decreasing particular solutions) 
and the continuous spectrum (solutions of the integral type) and that one 
particular solution disappears depending on the anomaly parameter. 

It has been established that the zeros of the dispersion functions are 
necessary for the analytical solution of the problem and, in particular, for 
deriving the electric field and the distribution function of electrons in the 
explicit form in the half-space. 

The analysis performed in this work has demonstrated that, in the case 
of the anomalous skin effect, the electric field in the skin effect problem is 
predominantly determined by the discrete spectrum. The real part of the 



electric field which corresponds to the continuous spectrum is eight orders 
of magnitude smaller than the real part of the electric field corresponding 
to the discrete spectrum. 

The imaginary part of the electric field at any anomaly parameters is 
considerably bigger than the real part. For £ = 10~^ the amplitude of 
imaginary part is four orders of magnitude bigger than the amplitude of 
real part. 
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Fig. 1 a. Domains in the {a, fl) plane. 




Fig. 1 b. Domains in the (uijUi) plane. 
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Fig. 2b. The imaginary part of the distribution function 
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Fig. 3b. The real part of the electric field for continuous spectra. 
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Fig. 3c. The imaginary part of the electric field for the discrete spectra. 
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Fig. 3d. The imaginary part of the electric field continuous spectra. 




Fig 4. The modulus of the electric field 
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